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Abstract. 

We present a numerical simulation study of the exchange bias (EB) effect in 
nanoparticles with core/shell structure aimed to unveil the microscopic origin of some 
of the experimental phenomenology associated to this effect. In particular, we have 
focused our study on the particle size and field cooling dependence of the hysteresis loop 
shifts. To this end, hysteresis loops after a field cooling process have been computed 
by means of Monte Carlo simulations based on a model that takes into account the 
peculiar properties of the core, shell and interfacial regions of the particle and the EB 
and coercive fields have been extracted from them. The results show that, as a general 
trend, the EB field Heb decreases with increasing particle size, in agreement with some 
experimental observations. However, closer inspection reveals notable oscillations of 
Jieb as a function of the particle radius which we show to be closely related to the 
net magnetization established after field cooling at the interfacial shell spins. For a 
particle with ferromagnetic interface coupling, we show that the magnitude and sign 
of Iieb can be varied with the magnetic field applied during the cooling process. 
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1. Introduction 

Magnetic nanoparticles have been the object of increasing research activity since the 
1950's. Their interest range covers from permanent magnetism, magnetic recording and 
magnetotransport [I], to macroscopic quantum phenomena or biomedical applications 
[2]. With the appearance of new preparation techniques, magnetic particles with sizes 
below the 10 nm range with narrow size distributions and different compositions can 
be prepared thorough a variety of procedures, either embedded in a conducting or 
insulating matrix or diluted in organic solvents. When entering the nanometer range, 
magnetic nanoparticles display interesting new phenomena that can be associated either 
to intrinsic characteristics of the individual particle (such as finite-size and surface effects 
[3]) or to collective properties of the nanoparticle ensamble, such as dipolar interactions 
[4] and exchange mediated coupling between the particles. 

In the last years, particular interest has been given to particles with core/shell 
structure of different compositions in which usually the core is ferromagnetic (FM) 
and the shell is made by the parent oxide compound [5], usually an antiferromagnet 
(AFM), although several studies have recently reported also inverse structures [6J. It 
is well known that proximity of a FM to an AFM induces exchange anisotropy at 
the interface between the two phases producing a shift of the hysteresis loop after 
colling in the presence of a magnetic field that has been termed exchange bias (EB). 
The origin of this effect has been intensively studied in thin film systems during the 
last decades and progress has been made in its understanding [7]. Although first 
reported for core/shell particles long time ago [8], there has been a renewed interest 
in understanding the microscopic mechanism underlying the EB effect in these systems 
after the observation that the EB effect can be used to beat the superparamagnetic limit 
of magnetic recording media [9j. Apart from the shift of the hysteresis loops, core/shell 
nanoparticles present peculiar magnetic phenomenology associated to EB which is not 
completely understood pUj due to the intrinsic inhomogeneity in the main parameters 
characterizing the core/shell structures present in real samples, such as particle size, 
shell thickness, anisotropy constants and easy axis directions. Theoretical models and 
theories suited for the explanation of EB in thin films bilayers are not directly applicable 
to core/shell nanoparticles and, in particular, are not able to account for the high values 
of the exchange bias fields Heb reported for some compositions [HI [12]. This is because 
the peculiar arrangement of spins at the core/shell interface of spherical or ellipsoidal 
nanoparticles results in intrinsic desorder and finite-size effects that are not present in 
the case of planar geometries. 

Recently, through numerical simulations of a microscopic model of an individual 
nanoparticle, we have shown that the bias field h^B can be quantitatively related to the 
net local exchange fields of the interfacial shell spins that act on the particle core [T3l[T4] . 
Our approach has also been successful in tracing back the microscopic origin of the loop 
assymetries observed experimentally to different magnetization reversal mechanisms 
in both loops branches [HI \W\. In the present study, we focus our interest on the 
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dependence of the EB effect on the particle size for a given value of th shell thickness 
and on the effect of the cooling field on the magnitude and direction of the shift of the 
hysteresis loops, unveiling its microscopic origin. 

2. Model and simulation method 

The Monte Carlo simulations reported here are based on the same model and 
methodology used previously in our recent works to which the reader is referred for 
the details [13]. The magnetic ions of the nanoparticle are modelled by classical spins 
Si lying in the nodes of a cubic lattice and interacting through the following by the 
following Hamiltonian: 

JV 

H/k B = Yl ■ % - k o E(^) 2 - k s E^ 2 ) 2 - E 1 ■ 3 w 

igC ieSh i=l 

where in the fourth term h is the magnetic field applied along the easy-axis direction 
with module h = fiH/k B in temperature units (/i is the magnetic moment of the spin). 
Spins in the different regions of the nanoparticle (core, shell and interface) have different 
magnetic properties. In the first term, the exchange constants are FM at the core 
(Jc = 10 K), AF at the shell (J$ = —5 K) and at the interface [formed by the spins in 
the core (shell) having neighbours in the shell (core)], FM or AF coupling Jj nt = ±0.5Jc 
will be considered. The second and third terms correspond to the uniaxial anisotropy 
energy, with k c = 1 and k s = 10 the anisotropy constants for core and shell spins 
respectively. We will consider particles with spherical shape and different total radius 
-Ribtai ranging from 6a to 20a and constant shell thickness Rsh = 3a (a is the lattice 
constant). The field cooling procedure previous to the simulation of the hysteresis loops 
has been performed from a temperature (T = 20 K) higher than the FM core ordering 
temperature down to T = 0.1 in steps 5T = 0.1 K and in the presence of a magnetic 
field hpc- The hysteresis loops have been computed cycling the magnetic field between 
h = ±h FC in steps 5h = 0.05 K, and the magnetization components averaged at every 
field value during 100 MC steps, discarding the initial 100 MCS after every field step. 

3. Results 

3.1. Particle size dependence 

The first objective of our work has been to study the dependence of the EB 
phenomenology on the particle size. In principle, the dependence of the exchange bias 
field /ieb should be similar to that observed in thin film systems on the thickness of the 
FM layer. If this was the case, one should expect an increase of /ieb when reducing 
the particle size(a e b ~ 1/Rc)- Some experiments have reported this trend in particle 
systems with different compositions [El ttH US] . However, other studies [2Ql EJJ reported 
the contrary variation in a certain range of particle sizes and some have even observed 
and argumented the disappearance of EB below a critical particle size [22j [23]. In 



Exchange bias in core/shell magnetic nanoparticles 4 



0.6 



^ 0.4 



0.2 
0.0 

-0.1 

a 

^ N -0.2 

-0.3 

-0.4 ( 

T(K) 

Figure 1. (Colour online) Thermal dependence of the magnetization on cooling in 
the presence of a magnetic field applied along the z axis hpc — 4 K for particles with 
Rsh — 3a and different radii R = 7, 9, 10, 12, 17, 20a (from bottom-up). The value of 
the shell and core anisotropy constants are fcs = 10 and fee = 1, respectively. The 
values of the exchange couplings are: Jc — 10 K, Js = — 0.5Jc, and Ji n t = — 0.5Jc- 
Panel (a) shows the total magnetization along the field direction, while panel (b) shows 
the contribution corresponding only to the interfacial spins at the shell for the same 
selected radii (same symbols as in panel a). 

order to clarify the origin of these discrepancies, we have first studied the magnetic 
configuration attained when cooling nanoparticles with different total radii R and the 
same shell thickness R$h = 3 in the presence of a constant field. The results of the 
cooling process are displayed in the figure CD, where we show the thermal dependence of 
the total magnetization (panel a) and also the contribution of the interfacial shell spins 
Mfj^ (panel b) for some representative particle sizes. The total magnetization increases 
with decreasing temperature, reflecting the progressive alignement of the FM core spins 
towards the field direction (see figure Uk). The values of magnetization attained at the 
lowest temperature increase with the particle size simply because of the increasing ratio 
of core spins with respect to the total number of spins in the particle. In contrast, 
the low temperature net magnetization of the AFM shell spins at the interfacial region 
varies in a nonmonotonous way with the particle size (see figure [lb)- Small changes in 
the particle core radius may induce very different geometric arrangement of spins at the 
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Figure 2. (Colour online) Dependence of the coercive field h c (panel a) and exchange 
bias field h c b (panel b) on the total particle radius i?Totai in lattice parameter units) 
for a shell thickness i?sh = 3a. The insets show a detail of the main panels around the 
region i?Totai = 11 - 13a. 



interface due to the particular intersections of the spherical shape with the lattice sites. 
As a consequence, the magnetic configuration of interfacial shell spins after field cooling, 
induced below the Neel temperatute of the AFM, is different for different particle sizes. 

In order to see the influence of these different magnetic states established after field 
coling for different particle sizes on the EB effect, we have also simulated hysteresis loops 
starting from the FC configurations for a range of particle size covering R = 6 — 20a 
and computed the values of the coercive field h c = (h^ — h^)/2 and the shift of the 
hysteresis loops h e \> = (^c + ^c)/^ from the values of the coercive field at the decreasing 
h~Q and increasing field Kq branches. The dependence of these two quantities on particle 
radius are shown in figure El We observe that, when decreasing the particle size, there is 
a trend for decreasing he which can be attributed to the higher proportion of interfacial 
core spins that have to be reversed along the hysteresis loop with increasing particle 
size. This observation, is in agreement with our previous finding [13] that the coercive 
fields after a FC process are higher than those obtained for zero field cooled hysteresis 
loops. For the EB field h e b, we observe that, although clear oscillations of this quantity 
with the particle size are seen, there is a trend to decrease as particle size increases. 
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The oscillations are in clear correspondence to the ones observed in M int in figure Q] 
and this is a demonstration that the microscopic origin of the loops shifts is related 
to the net magnetization component of the shell interfacial spins, since the local fields 
felt by the core spins are proportional to it, as we showed in our previous work [1 3j . 
Therefore, depending on the particular lattice geometry and the configuration of nearest 
neighbours at the interfacial spins, very different values of /i e b can be obtained, even for 
particles having close values of the radius. In order to demonstrate this, we computed 
additional hysteresis loops for some selected particle radii between 11a and 13a. The 
corresponding results for he and h e b are displayed in the insets of Fig. [2J, where the 
above mentioned fluctuations in both quantities can be clearly seen even for very close 
values of the particle radii. 

The reduction of the loop shift with particle size is in agreement with experimental 
observations reporting a similar trend [El [HI ES] . However, the opposite trend reported 
in some cases [20], [2Tj could be adscribed, within the scope of our model, to the fact that 
the limited range of particle sizes present in the samples of these experiments lain in 
between one of the oscillation periods found in our simulations giving rise to an increase 
in h eh instead of the expected decrease. However, it is possible that the existence of 
interparticle interactions and distribution of anisotropy easy-axis in real samples could 
also influence the dependence of h e \> on the particle size. 

3.2. Cooling field dependence 

It is well known from the experimental phenomenology on thin films formed by 
FM/AFM coupled bilayers, that the value and sign of the loop shifts may depend not 
only on the thickness of both layers but also on the nature (FM or AFM) of the coupling 
at the interface between them and on the value of the magnetic field applied during the 
FC process. For these kind of structures, it has even been demonstrated that loop shifts 
along the direction contrary to the cooling field (positive EB) can be obtained for high 
enough values of the cooling fields [24|. 

Fewer studies about the h-pc dependence of EB have been reported for core/shell 
nanoparticle systems, and there is no clear-cut interpretation of the results. While 
a study on Co/CoO nanoparticles [25] found an increase of h e h for h^c values up to 
5 T, in another study of Fe/FeO nanoparticles [26l [27] the authors observed a non- 
monotonic dependence with a maximum h e b at h^c = 5 T and a decrease above this 
value. In order to clarify these observations, we have performed simulations of hysteresis 
loops with different values of the cooling field h^c for a particle of radius i?Totai = 12a 
and FM interfacial coupling J\ nt = +0.5 Jo The results are presented in figure [3h for 
some representative values of /ipc — 4, 15, 20, 30 K, only the central part of the loops 
is displayed. In all cases, the loops are closed for values of \h\ lower than 4 K, but 
the shape of the loops becomes more symmetric with increasing hpc as indicated by 
the increased squaredness and higher values of the decreasing field branch remanent 
magnetization. However, as it is already aparent by direct inspection, the shift of the 
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Figure 3. (Colour online) Hysteresis loops for a particle of radius Rrotai = 12 and 
ferromagnetic coupling at the interface of value Jj„t = +0.5Jc obtained after field 
cooling from a high temperature in different magnetic fields Kfc = 4 K (red circles), 
15 K (green squares), 20 K (blue diamonds), 30 K (yellow triangles). 



loops (see figure [4b, where the dependence of on h FC is shown), which for low hpc is 
in the negative direction, becomes positive for h-pc > 25 K. In order to understand these 
observations, it is conveninent to look at the partiular contribution of the shell spins at 
the interface to the hysteresis loops, which are displayed in figure [3b, that present a clear 
dependence on /ipc- First of all, let us notice that the magnetization values attained by 
the interfacial shell spins after the FC process previous to the start of the hysteresis loop 
(Mj S n ^) increase with increasing /ipc, as signaled by the increasing magnetization values 
of the descending loop branch and by direct comparison of the values at the end of 
the FC process. Therefore, the effect of an increasing hpc is to reverse progressively the 
interfacial shell spins having one or two neighbours in the core into the field direction. 
The interfacial shell spins are the ones that generate the local exchange fields felt by 
the core spins. When cooling at low h^c, this exchange field points along a direction 
contrary to the magnetic field and produces a negative loop shift. Increasing h FC values, 
reduce the value of the local exchange fields leading to a reduction of the loop shift and 
ultimately inducing exchange fields that point along the magnetic field direction and 
loop shifts into the positive direction (see figure 0b). 

It should also be noticed that the decreasing and increasing field branches of the 
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Figure 4. (Colour online) Dependence of the coercive field h c (panel a) and 
exchange bias field h e b (panel b) on the magnetic field applied during the field cooling 
process previous to the calculation of the hysteresis loops for a particle with radius 
Rrotai = 12a. 

loops for Mf^ t in figure [3]b are quite different, with remanent magnetizations in the 
increasing branch lower (and even positive for the case h-pc = 4 K) than for the 
decreasing one. This difference reflects the fact that a fraction of the shell interfacial 
spins remain pinned during the field reversal, a phenomenon that is at the microscopic 
origin of the EB effect [T3j . As the cooling field increases, the difference between the 
remanent magnetizations decreases in accordance with the reduction of h e \,. Finally, the 
increase in coercivity with h-pc reported in figure Sh, is caused by the higher fraction 
of spins at the interface that become aligned into the core magnetization direction and 
that have to be reversed by the magnetic field dragged by the core spins. The maximum 
in the vs h-pQ curve reported in some works [26l IZZJ cannot be reproduced within 
the scope of our model. We think that, in order to explain this observation, additional 
features present in real samples, such as the policrystallinity of the AFM shell and 
surface disorder, should be incorporated into our model of nanoparticle, a work which 
is in progress at present. 
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4. Conclusion 

We have reported the results of simulations of the hysteretic properties of a model of 
a core/shell nanoparticle that have successfully explained the origin of some peculiar 
phenomenology associated to the EB effect. In particular, our results have linked 
the decrease of the loop shifts h e b with particle size to the microscopic magnetic 
configurations established at the particle interface after field cooling, which are 
responsible also for the oscillatory behavior of h e h with particle size found in the 
simulation results. We have also shown that, by varying the magnetic field applied 
during the FC process h eh , it is possible to modify the magnetic state of the spins at 
the interfacial region causing a corresponding decrease of /i e b and even loop shifts in 
the ositive field direction. We hope that further refinement of the model on which the 
simulations are based, by incorporating some source surface disorder will be able to 
account for training and time dependendent effects also observed experimentally. 

Acknowledgment s 

We acknowledge CESCA and CEPBA under coordination of C 4 for computer facilities. 
This work has been supported by the Spanish MEyC through the MAT2006-03999, 
NAN2004-08805-CO4-01/02 and Consolider-Ingenio 2010 CSD2006-00012 projects, and 
the Generalitat de Catalunya through the 2005SGR00969 DURSI project. 

References 

[1] Bader S D Rev. Mod. Phys. 2006 78 1. 

[2] Tartaj P Current Nanoscience 2006 2 43. 

[3] Batlle X and Labarta A J. Phys. D: Appl. Phys. 2002 35 R15. 

[4] Majetich S A and Sachan M J. Phys. D: Appl. Phys. 2006 39 (2006) R407. 

[5] Nogues J, Sort J, Langlais V, Skumryev V, Surinach S, Munoz J S and Baro M D Phys. Rep. 2005 
422 65. 

[6] Salazar-Alvarez G, Sort J, Surinach S, Baro M D and Nogues J J. Am. Chem. Soc. 2007 129 9102. 

[7] Nogues J and Schuller I K J. Magn. Magn. Mater. 1999 192 203. 

[8] Meiklejohn W H and Bean C P Phys. Rev. 1957 105 904 

[9] Skumryev V, Stoyanov S, Zhang Y, Hadjipanayis G, Givord D and Nogues J Nature 2003 423 850. 

[10] Iglesias O, Labarta A and Batlle X J. Nanosc. Nanotech. (in press, 2007) (Preprint 

|cond-mat/0607716| . 

[11] Kiwi M J. Magn. Magn. Mat. 2001 234 584. 

[12] Stamps R L J. Phys. D: Appl. Phys. 2000 33 R247. 

[13] Iglesias O, Batlle X and Labarta A Phys. Rev. B 2005 72 212401. 

[14] Iglesias O, Labarta A and Batlle X J. Phys.: Condens. Matter 2007 19 406232. 

[15] Iglesias O and Labarta A Physica B 2006 372 247. 

[16] Iglesias O, Batlle X and Labarta A ./. Magn. Magn. Mater. 2007 316 140. 

[17] Gangopadhyay S et al. J. Appl. Phys. 1993 73 6964. 

[18] Peng D L et al. Phys. Rev. B 2000 61 3103. 

[19] Del Bianco L et al. Phys. Rev. B 2002 66 174418. 

[20] Vazquez M et al. Physica B 2004 354 71. 



Exchange bias in core/shell magnetic nanoparticles 



10 



[21] Mumtaz A, Maaz K, Janjua B, Hasanain S K and Bertino M F J. Magn. Magn. Mat. 2007 313 
266. 

[22] Dobrynin AN et al. Appl. Phys. Lett. 2005 87 012501. 

[23] Martmez-Boubeta C et al. 2006 Phys. Rev. B 74 054430. 

[24] Nogues J, Leighton C and Schuller I K 2000 Phys. Rev. B 61 1315. 

[25] Zhou S M, Imhoff D, You-Zhang K and Leprince-Wang Y 2005 Appl. Phys. A 81 115. 

[26] Del Bianco L, Fiorani D, Testa A M, Bonetti E and Signorini L 2004 Phys. Rev. B 70 052401. 

[27] Fiorani D, Del Bianco L and Testa A M 2006 J. Magn. Magn. Mater. 300 179. 



